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1. Introduction 

The study of an open quantum system aims at understanding the underlying physical pro- 
cesses between a quantum system and the environment and at giving a correct account 
of the system noise due to the environment. Conventional theories usually assume a weak 
system-environment interaction. This assumption is no longer valid for a controlled quantum 
system, where the control, as a part of the environment, strongly interacts with the system. 
Quantum interference between the system and the control and decoherence environment can 
modify the noise in the quantum system. Given the recent advances in the experimental 
ultrafast optical techniques, an accurate theory that captures such an quantum interference 
effect is needed in the regime T2, where to and T2 are the operation and decoherence 
times, respectively. Also, the stringent error threshold requirement (e ~ 10~^) demanded 
by fault tolerant quantum computation [1] is another motivation for a high precision theory 
for open quantum systems. 

The master equation (ME) approach has been successful for Markovian open quantum 
systems [2], where the dynamics of the system is history independent. For such a Markovian 
system under a coherent control, the corresponding MEs are equivalent to the optical Bloch 
equations [3], where the control and dissipation terms are additive in the equations. How- 
ever, for systems with a structured environment [4], e.g. photonic crystal materials [5] and 
nanocavity systems [6, 7], the dynamics are non-Markovian and the optical Bloch theory 
is inapplicable. In fact, because of the memory effect of these systems, interferences exist 
between the control and the environment and the validity of the additive form of the ME is 
questionable. Therefore, the challenge remains open in non-Markovian systems. 

A precise diagrammatic technique has been formulated recently to solve the dynamics 
of a two level system (TLS) under a time-dependent coherent photon control through the 
multimode Jaynes-Cummings (JC) interaction in the regime to ^ ^2 [8]. It takes into full 
account of the quantum interference among the TLS, the control and the electromagnetic 
vacuum and is applicable for both Markovian and non-Markovian situations. The idea 
is to construct the time evolution of the TLS in terms of a time loop propagation, in 
the same spirit as the Keldysh non-equilibrium Green's function technique [9]. Under this 
construction, the underlying physical processes become transparent: dissipations originate 
from the contraction between photons at different times; coherent Rabi motions are due to 
propagators dressed by the control photons; interferences between these two processes result 
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in the control noise. For to <^ T2, only diagrams with a small number of contractions are 
important and they can be summed exactly, allowing a precise calculation of the controlled 
dynamics of the TLS in this time domain within a given error bound. Based on a simple non- 
Markovian model, it has been shown that this diagrammatic formalism captured quantum 
correlations of the entire system that was missed in conventional ME approximations even 
in the small time limit [8]. The error of the ME approximations is comparable to that of 
the classical approximation. 

This diagrammatic approach is not restricted to the case where the photon control is 
a Glauber state. The control can formally be extended to an arbitrary photon state. To 
illustrate this, we extend the theory to the cases of a multimode squeezed coherent state 
and a multimode photon number state. Comparisons with the single mode results indicate 
the accuracy of these diagrammatic solutions. 

The paper is organized as follows. In Section 2, we review the diagrammatic formalism 
for a TLS under a coherent control and provide rules for the diagrammatic construction. 
Section 3 performs an explicit comparison with the standard ME approaches. We also 
demonstrate the ability of the diagrammatic approach to systematically improve the accu- 
racy by including higher order diagrams. Section 4 extends the discussion to non-Glauber 
photon controls. We conclude in Section 5. 

2. TLS under coherent control 

The central idea of the diagrammatic technique is first to express the time evolution of the 
composite system in terms of the spin and photon operators, then to evaluate the transition 
matrix in terms of photon operators, which are finally treated diagrammatically using the 
Wick's expansion. The problem of a time-dependent controlled TLS is governed by the 
multimode JC Hamiltonian under a multimode coherent state \a) = \aki, a^a, • • •)• Consider 
first for simplicity, an initial wavefunction of the product form |\1'(0)) = Cs\s)] \a.), where 
s = ± gives the spin state. The system evolves as |\l'(t)) = f/(t)|\I^(0)), where the evolution 
operator in the interaction picture is given by: 




(1) 
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and the multimode JC interaction is: 



V{ti) = a+Ai + a^Al (2) 
with Ai = J^dkakc'^'^'-''''^''. (3) 

k 

with uq the energy sphtting of the TLS, Uk the photon frequency and gk the TLS-photon 
couphng strength. The interaction on the initial coherent photon state can be viewed as the 
action of a control pulse, consisting of a modulation of the resonant frequency mode with 
the coherent state amplitude satisfying: 

{9ke-'-^W + c.c) = ^e-^(-o*+^) + c.c, (4) 

where is the eigenvalue of on \ct), Q(t) and are the envelope function and phase of 
the pulse, respectively. Note that we have made the rotating wave approximation which is 
valid when u^t g\a\t 0(1) [8]. The physical quantity of interest is the transformation 
matrix: 

Psf,sy,sAi,cx) = H{s'\UHt)\s'f){sf\U{t)\s)\(x) 

from which we can obtain the reduced density matrix of the TLS at that: 

Psj,s'j{t) = Ycscl,ps^,sy,s,s'it,cx). (6) 

s,s' 

It is straightforward to generalize the above formulae to an initial state which is a linear 
combination of the product states. 

In contrast to the ME methodology which traces out the photonic environment in the 
beginning, we evaluate the TLS transition by the spin operators in the transformation matrix 
first, leaving a sequence of the corresponding photon components. An example is given by: 

P++;++{t) 



n,n'=0 



JO 



X (a| [AvAl,...A2n'-lAl^}j (^A2„4n-l-^24) l«) , 

(7) 
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where Jq D'^t = J^dtn... J^^ dt2 J^^ dti. In the following, we will focus on this expression for 
illustration. Other transformation matrices can be obtained similarly (see [8]). Therefore, 
the task is now reduced to compute expectation values of the photon operators for the initial 
multimode coherent state. 

Eq. (7) is exact up to now and all the quantum correlations of the whole system remain 
intact. The underlying physical processes can be easily identified using a diagrammatic 
representation as in Fig. 1. We shall see in the next subsection how different roles of 
the photons affect the interaction dynamics, and how this infinite series expression can be 
largely simplified and summed to a desirable accuracy using t/T2 as a small parameter. No 
stochastic assumption has been made. We emphasize that this approach is fundamentally 
different from the ME method which usually assumes a stationary environment and thus 
does not account for the interference effect between the control and the environment of a 
quantum system. 

A. Diagrammatic solution 

It is instructive to describe Eq. (7) as time loop integrals of photon operators propagating 
forward and backward between time and t, and depict Fig. 1 in a bubble form. Diagram- 
matic series is constructed using the Wick's theorem [10], which expands the product of 
photon operators in terms of the normal ordered form and possible contractions between 
two photon operators. By doing so, the expectation value in Eq. (7) can be computed easily 
by the replacement — )■ ak and a\ — )■ for normal ordered photons. Physically, the 
normal ordered photons coming from the multimode coherent state are responsible for the 
coherent Rabi motion of the TLS; while the contraction between photons, being independent 
of the control, corresponds to vacuum relaxation of the TLS. 

The diagrammatic structure consists of two ingredients: the contraction between photons 
and dressed lines due to the coherent control. The contraction function is given by: 

{AA]) = K{t, - t,) = (7^e^(-o-^)(*-*.), (8) 

k 

and is depicted in Fig. 2a. Here, ti and tj can be on the same or opposite time lines. Note that 
the contraction only depends on the photon density of state (DOS) p{uj) = J2k9k^i^ — <^k)- 
One can show that each contraction line ^ j J dtidtjK{ti — tj) ~ 0[{t /T2y']., where 1 < 
7 < 2 depends on the DOS [8]. Therefore, the Wick's expansion in terms of contraction is 
automatically a perturbative series of the small parameter {t/T2)^ . On the other hand, the 
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dressed lines lead to Rabi oscillations of the TLS and are given by [8]: 

D^{t,t') = {±^)ie^"^) sin ©(^ _ t'), (9) 

corresponding to the double and triple line depicted in Fig. 2b. A dressed line segment 
between the same (opposite) spin state corresponds to a double (triple) line. Here, Do{t,t') 
picks the factor +i {—i), when the triple line is on the upper (lower) time line; and gains 
a phase e^'^ (e"*"^), if it goes from — to + in the clockwise (anticlockwise) sense. A(t) = 
dt'Vtit') gives the area of the coherent pulse at time t. 

The diagrammatic procedure can be summarized as follows: (i) denote the initial and 
final state of the transformation matrix on a loop diagram, where the lower and upper arcs 
represent forward and backward time evolutions; (ii) draw all possible contractions according 
to Fig. 2a; and (iii) dress the bare lines according to the spin flips by the photons. Each 
diagram can be computed straightforwardly by integrating over all time variables at the 
vertices of each contraction line. The vertex picks up a factor of i {—i) when it is on the 
upper (lower) time line. 

The diagrammatic structure allows an easy identification of the underlying coherent and 
dissipative processes. For instance. Fig. 2c shows the diagram with no contraction for 
p_i_+._i_+(t). Following the above description, this diagram simply gives the expected coherent 
Rabi solution p^+.^+(t) = cos^ [^(t)/2]. On the other hand, in the absence of control, the 
vacuum relaxation is represented by a series of undressed diagrams that involves contractions 
only, as can be seen in Fig. 2d. In the regime t <^ T2, the first two diagrams in Fig. 2d 
results in a decoherence of ~ 0[{t/T2)^]. 

The leading order contribution to the control noise problem is represented by the three 
diagrams in Fig. 2e. Interference between the control (dressed line) and decoherence (con- 
traction line) is evident. We stress that these diagrams are first order in decoherence, but 
infinite order in the coherent interaction. This explains why the diagrammatic approach is a 
suitable candidate to our control noise problem, in contrast to the ME method that requires 
a weak system-environment coupling. 
B. Application to a general DOS 

The field theoretic approach described above is applicable for an arbitrary photon DOS and 
thus valid for both the Markovian and non-Markovian systems. For a broadband DOS, the 
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contraction function becomes a delta function K{ti — tj) ~ S{ti — tj), meaning that the two 
ends of a contraction hne are squeezed to the same point. In this hmit, the correlation time 
Tc is zero and the decoherence is linear in t/T2 (i.e. 7 = 1), corresponding to a Markovian 
exponential decay. The dynamics of the controlled TLS can be alternatively obtained by the 
optical Bloch equation and the diagrammatic solution shows quantitative agreement with it 
in the regime t <^ T2 [8]. 

Extension to the non-Markovian regime can be performed simply by modifying the photon 
DOS and thus the contraction function. The correlation time Tc is characterized by the 
bandwidth of the DOS. The contraction function is no longer a delta function, leading to 
non-exponential decoherence ~ 0[(t/T2)'^], where 1 < 7 < 2. Furthermore, there is a 
transition from the non-exponential to exponential decay as t becomes comparable to r^,. 
All these features remain in the presence of a coherent control [8]. In the following section, 
we will investigate the accuracy of the diagrammatic scheme in comparison with existing 
ME methods. 

3. Comparison with ME approaches 

Conventional master equation (ME) assumes an additive form [11]: 

j^Psit) = -i[H,{t),pM + f^dt'L{t-t')p,{t'), (10) 

where Hc{t) describes a control Hamiltonian that only acts on the system and L{t — t') is 
a suitably chosen superoperator that only accounts for the dissipative environment. It is 
assumed that these two terms are independent on each other. In this section, we examine 
the validity of this kind of ME against the field theoretic solution in Section 2. 
A. Second order TLS- environment interaction 

We consider two commonly adopted ME approaches based on the projective operator tech- 
niques, namely the Nakajima-Zwanzig (NZ) and the time-convolutionless (TCL) methods 
[12, 13]. We now show that they belong to the additive ME in Eq. (10). The second order 
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NZ ME is given by: 




+ / dt'TiR [Vit), Ttr [V(t'), pf (f) ® Pi?] ® Pr] 



tdt'Tm [Vit),[Vit'),p^'it')0pn\], 



(11) 



Jo 



where ps{t) is the reduced density matrix of the TLS, p/j = |q:)(q:| denotes the initial 
multimode photon density matrix and has no time dependence. The non-vanishing first and 
second terms are due to fact that TT^[V{t) (g) p^] ^ 0. Using the multimode JC interaction 
in Eq. (1) and the relation in Eq. (4), we trace out the photon reservoir and arrive at the 



where K{t — t') is just the contraction function in Eq. (8). The second order TCL ME can 
be obtained simply by the replacement p^^^{t') — )■ p^^^{t) in the integrand. 

It is clear from Eq. (12) that the NZ ME and TCL ME belong to the class of additive 
ME in Eq. (10), where the corresponding superoperators L(t — t') do not depend on the 
control, and thus both of them have neglected the quantum interference effect between the 
control and dissipation. This can be explicitly illustrated by using a simple non-Markovian 
system, the single mode JC model, where an exact solution is available. This situation 
corresponds to a constant contraction function K{t — t') = g^, infinite correlation time 
Tc = oo, 7 = 2 and T2 = V^/g- We consider an initially excited TLS under a single 
mode coherent drive with n = lOOvr^, so that the system undergoes a Att rotation when 
gt = 0.2. Fig. 3 compares the difference between different approaches and the exact solution 
for P++(t) = (t). Here, the diagrammatic solution, based on the three diagrams in 

Fig. 2e, shows an excellent agreement with the exact solution in the small gt regime. The 
error of these three diagrams is 0[{gt)'^]. On the contrary, these two ME methods, expected 



NZ ME: 



d 
di 





(12) 
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to be the same and accurate to the second order of gt [12, 13], show errors of the same 
order of magnitude. The error of the NZ solution is very close to that of the classical Rabi 
solution (P++(t) = cos^((7v^^)) which entirely neglects any decoherence effect; while the 
error of TCL result shows a phase difference from the NZ and classical solution because of 
the approximation p^^^{t') — > p^'-^^{t) in the integrand of the ME. 

The problem of these ME approximations stem from the assumptions of a stationary 
environment, which is inappropriate for a quantum system under a general control. Quantum 
feedback from the photon environment to the TLS exists and originates from the quantum 
interference between the control and the environment. In fact, the difference between the 
field theoretic solution and the additive ME method can be viewed as the difference between 
the quantum and semiclassical correlations. This difference is not limited to the single mode 
test above and occurs for a general non-Markovian multimode system. 
B. Higher order terms 

The precision of the diagrammatic solution can be further improved by including higher 
order diagrams. In general, the inclusion of diagrams with n contractions would result 
in an error of 0[(t/T2) "'^^'^]. Fig. 4 provides the next order diagrams that possess two 
contractions for p^^.^^^{t). They can be evaluated by using the same procedure described in 
Section 2 A. As an example. Fig. 5 plots the absolute difference of P++(t) between different 
methods and the exact solution in log scale under the same physical situation as in Fig. 3. 
It is apparent that the diagrammatic solution has a well-controlled error bound in the small 
time domain. 

4. Control by photon states other than coherent state 

The field theoretic technique can be formulated similarly when the TLS is initially in the 
presence of a quantum photon state other than the coherent state. The Wick's expansion 
is still valid, though the diagrammatic rules have to be modified. Here, we extend this for- 
malism to the multimode squeezed coherent state and the multimode photon number state. 
Rather than exploring all possible scenarios, the purpose of this section is to demonstrate 
how diagrammatic solutions can be constructed with different initial photon states. 
A. Squeezed coherent state 

A multimode squeezed coherent state takes the form |q:,^) = D{cx)S{S,) |0), where D{a) is 
the multimode displacement operator and we take the multimode squeezing operator to be 
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S = rifc ^^P yC'^k+k'^k-k ~ [14]. k governs on the squeezing mechanism. The 

problem is to evaluate Eq. (7) using \cx.,^) instead of \a). 

Making use of the squeezing transformation, we can rewrite Eq. (7) by the replacement 
Ai^ Ai = S^AiS, a] A] = S^A]S and D{a) D{a) = S^D{a)S, where each photon 
operator transforms according to 

5*^0^5 = flfc cosh r — a,2fc-fc^*^ ^^^^ 

S^a^S = a|, cosh r — a2A;-fce~*^ sinh r, (13) 

and the squeezing parameters satisfy ^ = re'^. Now, we can proceed the Wick's expansion 
as in Section 2 A. The dressed lines are unaffected as the normal ordered photon operators 
are just c-numbers. As such, the zeroth order diagram remains the same. However, four 
types of contractions are now possible due to the squeezing: 

{AiA]) = coshV^f^^giK— 

k 

(AlAj) = sinhV^(72g-i(^o--.)(t,-ti)^ 

k 

(AA,) = -— sinh(2r) 

k 

{AIA]) = {AA,y, (14) 

and they are depicted in Fig. 6. The above dependence of the contraction functions on the 
squeezing parameter 6 shows a phase sensitive decoherence of the TLS. 

In brief, the diagrammatic rules in Section 2 A still apply with the exception that the 
construction of contractions in Fig. 2a is replaced by that in Fig. 6. The leading order 
contribution to control noise now consists of twelve diagrams which are similar to the three 
shown in Fig. 2e. To demonstrate the accuracy of the diagrammatic solution, we compare 
it with an exact solution of a single mode JC system under a squeezed coherent state 
[15]. Fig. 7 again shows an excellent precision of the diagrammatic solution. Depending 
on the squeezing parameters, the phase sensitive decoherence of the TLS can be greater or 
smaller than that without squeezing. We note that since the decoherence time is modified 
to ~ l/(5fcoshr), the small time requirement becomes 5fcosh(r)t ^ 1 for the validity of the 
perturbative solution. 
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The above squeezed state analysis works similarly for the more general linearly trans- 
formed coherent state D{(x)T \0), where T is any unitary operator that transforms the 
photon operators linearly, T^^a^T = J2k' i^kk'Q'k' + Vkk'O^k^ ■ The diagrammatic solution 
shows that the coherent motion of the TLS remains the same, while the decoherence can be 
computed straightforwardly by using the modified contraction functions similar to Eq. (14). 

B. Number state 

The dynamics between a TLS and a multimode number state is more complicated. A single 
photon process can project the photon state to be orthogonal to the initial state. It is a 
daunting task to keep track of the time evolution of a number state with an arbitrary photon 
number distribution. With the help of the diagrammatic technique, we shall see below how it 
is possible to tackle a number state of the form |A^, {0}), where {0} is a multimode vacuum. 

The idea is again to make use of the Wick's expansion of the photon operators in Eq. (7), 
where the coherent state is replaced by |A^, {0}) here. We assume the A^-mode of the number 
state has zero detuning. The zeroth order term has all the photon operators being normal 
ordered, so that each photon operator can only annihilate or create photons in the A^-mode. 
Thus, the time integrals in Eq. (7) can be performed and the resultant expression is: 



which is already different from the classical Rabi solution. The coherent photon dressing 
described in Section 2 A no longer applies and one has to perform an infinite sum as in 



The contraction line remains as a good small parameter. The first contraction terms 
can be calculated as in the way to obtain Eq. (15). Fig. 8 shows the relaxation, defined 
by P++(t) — p^"^*^(*)^ of a TLS driven by different initial photon states. For simplicity, 
we consider the vacuum mode to be degenerate with the A^-mode, so that the contraction 
function is time independent. We define the driving term P^"™(t) = cos^ {gy/N + 1 t) and 
cos"^ {g\a\t) for the number and coherent states, respectively. In the case of a single mode 
number state |A^), no relaxation is seen as expected. However, when we consider a number 
state with one degenerate vacuum mode |A^, 0), dissipation is observed due to the interference 
between the control photon and the vacuum. The amount of decoherence increases linearly 



oo 

= ^-N{gtr+ (-1) 

n,n'=l 




(15) 



Eq. (15). 
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with the number of degenerate vacuum modes present. It is also interesting to observe that 
the decoherence of |A^, 0) is equal to that of a single mode coherent state |a) in the small 
time regime. 

The diagrammatic approach serves as a first step towards a precise treatment of this 
multimode number state problem. The generalization to deal with states having a general 
photon number distribution remains as a challenge. The ME approximation schemes are even 
worse at accounting for the decoherence effect of this problem compared to their efficiencies 
in coherent state case. Following the discussion in Section 3, one can show that for an initial 
photon state |A^, 0), the NZ ME leads to a non-dissipative oscillations of the TLS, where 
P++(t) can be negative even in the small time regime gt <^ 1. The TCL ME results in a 
Gaussian decay {P^^^{t) ^ [1 + exp {—2NgH'^)]/2, for ^ 1), where the Rabi oscillations 
are lost. Their failures again originate from the inability to capture the coupled dynamics 
between the control and the environment. 

5. Conclusion 

The diagrammatic formalism provides a new platform to understand and compute the inter- 
action dynamics among the quantum system, the control photon and the environment. Its 
high accuracy is attributed by the full account of the unitary evolution of the wavefunction 
of the whole system without any stochastic assumption. The revelation of underlying phys- 
ical processes by the diagrammatic structures enables a perturbative computation scheme 
for systems under fast photon control versus slow relaxation, with well-defined error bounds. 
It also serves as an extension of the standard field theoretic techniques in quantum electro- 
dynamics and many-body problems to systems far from equilibrium. 

The ME approach has been a good phenomenology to provide qualitative understandings 
of open quantum systems. Given the recent advances in high precision quantum technologies, 
there is a need to reexamine and improve the accuracies of the existing theories. The 
diagrammatic solution indicates the absence of quantum interference between the control 
and the environment in current ME approximations. This general effect is not restricted to a 
TLS interacting with a single mode photon state. We believe that the diagrammatic results 
can motivate improvements of conventional ME studies and stimulate the developments of 
other novel quantum theories. There is hope that more quantum mechanical effects of open 
quantum systems will be uncovered in the future. 
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Fig. 1. An illustration of the TLS-multiphoton process conditioned on a multimode coherent 
state. The ± sign labels the spin state of the TLS. Contractions among photons (dashed 
line) correspond to decoherence, while the normal-ordered control photons (wavy lines) is 
responsible for the Rabi motion of the TLS. diagram0.eps 
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Fig. 2. Diagrammatic representations of the Wick's expansion of the transformation matrix 
Psf,sy,s,s'{t, {<y}) for a coherent control, (a) The two possible contractions between photons 
at different times, (b) Two types of dressed line that represent a sum of even or odd numbers 
of control photons interacting with the TLS. (c) A dressed diagram without any contraction 
represents the classical Rabi solution of from time to t. (d) Vacuum relaxation 

is represented by undressed diagrams with contractions only, (e) The leading contribution 
of control noise to comes from three dressed diagrams with only one contraction. 
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Fig. 3. (color online) A comparison between the diagrammatic methods and the ME ap- 
proaches presented in [8], using a coherently driven single mode JC system with n = lOOvr^, 
corresponding to a Air rotation at gt = 0.2. The magnitudes of errors of the ME approaches 
and the classical Rabi solution are comparable even in the small time regime, compari- 
son!. eps 
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Fig. 4. Fifteen dressed diagrams with two contractions for under a coherent controL 
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Fig. 5. (color online) A plot of absolute errors of different methods under the same phys- 
ical situation as in Fig. 3. Diagrammatic solution with n contractions has an error of 
O [(s't)^'-'""*"^^] . comparison21egend.eps 
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Fig. 6. Eight possible contractions between photons for a squeezed vacuum, diagram3.eps 
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Fig. 7. P++{t) - Pf^^'^^^t) of an initially exited TLS under a single mode squeezed coherent 
state with n = lOOvr^ for various squeezing parameters, (p is the phase of a* and 6 is the 
phase of squeezing. The curves and symbols correspond to the diagrammatic and exact 
solutions, respectively, squeezedstate.eps 
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Fig. 8. Relaxation of an initially excited TLS under different initial photon states. 
N = lap = 103. The case with a non-dissipative single mode number state (|A^)) gains 
decoherence when a degenerate vacuum mode is added (|A^, 0)). In the small time region, 
the amount of dissipation is the same as that of a single mode coherent state num- 
berstate.eps 
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